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1. INTRODUCTION 

Renewable energy is becoming a popular source of clean energy these days [1]. Due to its cost 
efficiency, trivial pollution and continuous supply, they are becoming the latest trend [2]. Wind energy has 
become a primary non-conventional source of energy with increasing installed capacity and offshore 
generation [3]. In recent times, efforts have been made to exploit the major environmental benefits of wind 
energy reaching up to 5.5x10°-MW of world installed capacity by 2017 and increasing offshore installations. 

In today’s time, grid connectable systems have assumed importance as they directly export power 
for usage and void storage necessities [4]. In a grid connected, non-conventional system such as wind 
turbine, the input energy is intermittent and hence continuously changes with time. Correspondingly, 
the output power in such system changes with respect to time at every instant. In such a grid connected area, 
system instability, faults and transients may damage our turbines due to electromagnetic stresses [5]. 
Therefore, besides generating maximum power from renewable sources, efficient generation of energy is 
a prime requisite [6]. 

With the continuously varying wind speed, the output active power also changes. A slight variation 
in the wind speed may cause the system to become volatile to high forces and speed causing friction in 
mechanical parts. This demands power control and parameter adjustments which is done by using automatic 
controllers. Recent studies have proposed a wide range of control algorithms [7]. A PQ controller is used to 
improve system performance in renewable energy systems to bring output power to a nominal value [8]. 
However, implementation of such controllers 1s complex. An effective method for system control is fuzzy 
logic control instead of conventional PID controllers [7]. A combination of fuzzy inference system and neural 
networks has evolved as an optimal control method. However, the adjustment of fuzzy membership functions 
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through hit and trial is laborious. Therefore, ANFIS is being used for systematic construction of fuzzy rules 
and membership functions [9]. 

Besides active power, a system needs reactive power control. With the change in wind speed, 
the reactive power demand of the system also changes. And therefore, grid control is becoming a necessity 
with the continuous change in reactive power requirements of the system. An induction generator is excited 
by an externally connected capacitor bank. The capacitor size varies according to system needs and 
accordingly becomes bulky and costly and can serve a fixed reactive power demand. Therefore, now-a-days 
STATCOM is being used for variable reactive power compensation. Besides being capable of providing 
variable reactive power compensation, STATCOM is a cheap alternative to operator requirements and 
variable power needs [10, 11]. A combination of a static capacitor bank and STATCOM serves as an ideal 
choice for systems where sporadic behavior of input is expected [12]. 

Continuously varying inputs and surges may introduce transients in the system. It is important to 
carry out transient analysis to maintain overall security of the system and keep system in operation. Transient 
stability studies and system monitoring helps in extracting productive energy by limiting voltages, currents 
and speed. Therefore, in renewable energy systems like wind turbines, transient study assumes prime 
importance. Transient state analysis can be well performed with MATLAB/SIMULINK environment and 
MATLAB environment permits us to carry out varied analysis through simulation methods: 

-  Phasor method 
- Discrete method 
- Continuous method 

All of them have their respective applications and limitations thereby judging their usage in any 
model. A continuous simulation method is applied for an accurate solution in systems having continuous 
state. Whereas, a discrete model is used in a discrete model where a state changes at a particular time and 
remains for some time. Phasor method finds its usage in transient stability analysis for large power systems, 
motors and generators with long simulations times. It is used to study oscillations and to find simpler 
algebraic equations representing current and voltage phasors. Nonetheless, it can be applied to linear circuits 
to calculate phasor changes in voltage and current. It is a fast method however, can be applied only at one 
frequency [13]. 

In this paper, the electrical and mechanical performance of an IG wind farm has been improved 
by implementing an ANFIS based controller in MATLAB/SIMULINK environment. Since the scope of 
the paper is to conduct transient state study, phasor analysis has been appropriately used to achieve the same. 
Further, after obtaining results with PI and developed controller, the detailed system model is linearized and 
stability of system verified using pole-zero plot. 


2. SYSTEM DESCRIPTION AND ITS CONTROLLERS 

A 9-MW wind farm [13], consisting of three 3-MW generation pairs, with its control block has been 
depicted in Figures 1(a) and 1(b) respectively. The wind turbines are coupled with squirrel cage induction 
generators. The system is connected to 120 kV, 60Hz grid through a 25 km, 25 kV feeder. The stator winding 
is connected to the grid and a variable-pitch wind turbine drives the rotor. The pitch angle control limits 
the power flow when the wind speed exceeds the nominal value of 10m/sec. The speed ranges between 1 pu 
and 1.005 pu for unloaded and fully loaded condition respectively. The detailed system parameters have been 
given in Table 1. The system consists of four basic components i.e. wind turbine, induction generator, 
reactive power compensation block and control block, each of which has been described. The values of 
parameters for wind farm are given as shown in Table 1. 





(a) 
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Figure 1. (a) Diagram representing control block, (b) Grid connections with IG wind farm 


2.1. Wind turbine 

Wind turbines convert wind energy into mechanical power. Wind turns the turbine blades, 
connected to rotor, which spins the generator to generate electricity. The mechanical output power is given to 
the induction generator. It converts mechanical output power into electrical power. As discussed earlier, 
the mentioned system is a 9-MW wind farm consisting of three 3-MW generation sets. Each of the generator 
set comprises of two 1.5-MW turbines to export power to 120-kV grid. The mechanical output power for 
the system is given as [14]: 


P mech. ey yLareaP Vi 
Where, 
Mechanical Output power = P meh 
Performance coefficient = Poerf 
Air density = 
Pitch angle of blade =y 
Area swept by blade Z 
Wind speed =V% 
Tip speed ratio =p 


2.2. Induction generator 

An induction generator is an alternating current (AC) machine that operates on the principle of 
induction motors to produce power. The rotor runs faster than synchronous speed with the help of a prime 
mover. The wind farm uses Squirrel-cage induction generators. Squirrel-cage generators are most commonly 
used in induction generator applications. They consist of a squirrel-cage rotor which is a cylinder with 
embedded Aluminium and Copper conductors. 


2.3. Reactive power compensation 

Each pair of 1.5MW turbine demands 400-kvar of reactive power. During steady state, the capacitor 
banks are capable enough to support reactive power generation. However, during transient state, 
as the reactive power demand increases, it is compensated by STATCOM. STATCOM is an economic 
solution to present compensators, around 40% cheaper than its alternatives [6]. In absence of STATCOM, 
the voltage at bus drops causing overload of the wind turbine. It maintains a voltage of 25 kV at the 
connected bus during transient condition. 
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2.4. Control block 

With the change in system active power, the system demands power control and adjustments. This is 
achieved by the control block (refer Figure 1(a). The change in power Pn is fed to the a comparator and 
an error signal P, is generated by comparing it with a reference value P ef. The error signal is fed to 
the controller for system control. The controller is a conventional PI controller which is then replaced 
by ANFIS based controller to improve system performance. Both PI and ANFIS based controllers have been 
developed in the next sections to come. With the change in wind speed, the pitch angle change is a, for PI 
controller and a, for ANFIS based controller. 


3. DESIGN OF PI CONTROLLER 

A conventional PI controller is to be connected for system control. The model is subjected to linear 
analysis mechanism under control design. For this, the linearization tool is used. The methodology allows us 
to gather system information and parameters. At first, with linearization of model, we get system transfer 
function which is as follows: 
Let the controller structure be represented by T,. 


Where, (1-2) 
_ num 
c den 
num = a s nta a Pennin +atst+a?s® (1) 
den = p76 eso va sisesbastuse +bts1+b°s° (2) 


The value of all coefficients have been presented in Table 2. 
The controller structure is of the type kpro + wins [15]. 


T, can also be written as: 
Consider the plant transfer function to be T». 
The transfer function can also be written in exponential form as [16]: 


L 0 
T, = Mpeak€e J 
Therefore, the transfer function for open loop system becomes T,T,, . 


Let the phase margin for given system be @ and open loop crossover frequency to be. 
With the given phase margin and crossover frequency, open loop transfer function can be given as: 





(a [im ) (Mpeak ej®) — pi(-180+9) 

To evaluate kpro and kint , we solve above (3-4): 
-kpro OMpeax sin sin 0 + kintMpear COS cos 8 = w sin sin Ø (3) 
kpro OMyeax COS cos O + KineMpeagx Sin sin 0 = —w cos cos Ø (4) 


From above equations, we find the following (5-4): 


—cos(9-8) 


Kpro = coe (5) 
kine = SA (6) 
peak 


With the knowledge of Mpeag» œ, Ø, 0 , we can calculate the value for kpro and kint. The values of 
which have been given in Table 1. Whereas the values of coefficients are given in Table 2. The design of PI 
controller is essential for design of the ANFIS based controller. The design of the proposed ANFIS based 
controller is given in the next section. The various system parameters for model under consideration are 
as given Table 1 and 2. The values for coefficients in equations (1) and (2) are given in Table 2. 
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Table 1. Parameters of 9-MW wind farm 


Turbine Output Power 3 MW each 
Generator inductance 6.77 pu 
Turbine Base Wind Speed 10m/sec 
Koro 5 
Kint 25 
Maximum Pitch Angle 45 deg 
Fixed Capacitor reactive power 400Kvar 


Table 2. Table for values of coefficients 


NUM DEN 
Coefficient Value Coefficient Value 

ae? 81.61 b” 1 
ar 1.267e05 pee 1599 
a? 2.366e08 b?° 3.055e06 
a“ 2.442e11 b2* 3.238e09 
az! 2.501e14 pas 3.39e12 
a’ 1.887e17 pe? 2.618e15 
a 1.323e20 bp% 1.874e18 
a? 7.617e22 b2° 1.101e21 
ai’ 3.899e25 bt? 5.756e23 
gs 1.746e28 p18 2.62e26 
a 6.567e30 bY 1.008e29 
av 2.281e33 pe 3.543e31 
at’ 6.051e35 bi 9.631e33 
qi 1.582e38 bi4 2.534e36 
ar 2.594e40 pe 4.286e38 
at? 4.579e42 bY 7.56e40 
a? 2.852e44 pb 4.969e42 
aè 9.907e45 b1? 1.879e44 
a’ 2.365047 b? 4.976e45 
a® 3.547e48 b® 8.942e46 
a’ 4.069e49 b7 1.259e48 
af 3.273e50 b® 1.331e49 
ae 1.944e51 b” 1.12e50 
a’ 7.83e5 1 b* 7.165e50 
at 1.716e52 b? 3.465e51 
a? 1.438e52 b? 1.184e52 

b1 2.337e52 

b? 1.841e52 


4. DESIGN OF ANFIS BASED CONTROLLER 

Fuzzy logic control is gaining wide significance as a control technique. Fuzzy logic controllers are 
gaining importance in terms of improving turbine aerodynamics, improving efficiency and enhancing 
performance. However, it cannot automatically tune the set of rules. Whereas adaptive systems provide 
comprehensive training mechanism and can be applied to a wide range of practical systems. 
This combination of fuzzy inference and adaptive neural networks has emerged as a novel computational 
method against all of its counterparts [16]. A detailed study on ANFIS has been given in [17, 18] and [19]. 

The working of adaptive neural methodology is similar to neural networks. It provides information 
on a set of data by computing membership function parameters, mapped to a set of rules. It maps rules to 
output characteristics and output characteristics to single valued output or decision. The system is modelled 
by the fuzzy inference constructed by ANFIS with tuned parameters and training algorithms. Adaptive 
learning through ANFIS consists of four marked spheres of parameter tuning. 
- Parameter adjustment 
- Data Identification 
- Model validation 
- Testing data sets 

ANFIS editor GUI helps us to train and test Sugeno-type fuzzy inference system. The editor can be 
accessed through a MATLAB prompt anfiseditor. The subsequent tasks in ANFIS editor have been explained 
in [14]. To train the controller with ANFIS GUI, the data is loaded from the PI controller into the editor 
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interface. The training data set are collected on the basis of observations and target system. The output of PI 
is fed into the editor for training ANFIS. Flow chart representing GUI tasks as shown in Figure 2. 

A structure can be developed from the loaded data by varying the number of membership functions. 
The data is repetitively trained according to varied membership functions so as to get an accurate training 
result. After each training, simulations are run and corresponding data is observed. Figures 3-7 show 
the developed membership functions, structure and rules obtained for an accurate simulation. Fuzzy if-then 
rules are associated with membership function used during training .The fuzzy if-then rules for given 
simulation as shown in Table 3: 

- If G@nputl is inlmf1) then (output is out] mf1) (1) 
- If (input! is inlmf2) then (output is out] mf2) (1) 
- If G@nputl is inlmf3) then (output is out] mf3) (1) 

After subsequent training, the data is tested on data sets on which it was not trained. By validating 
the model, we check the general working capability of Fuzzy logic control on the given system. The scope of 
the paper is to implement the principle of balance. The control action is performed in order to maintain 
a nominal value of active power. A PI controlled system was able to maintain power control by changing 
the pitch angle. However, with the help of ANFIS based controller, a better control action can be performed 
by a smaller change in pitch angle thereby improving system performance. 


Table 3. Specifications for ANFIS 


Parameter Value 

No. of inputs 1 

No. of outputs 1 

No. of membership functions 3 
Type of membership function trimf 

No. of epoche 3 
Average testing Error 1.2842 











Figure 2. Flow chart representing GUI tasks 
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Figure 5. Membership functions for trained data Figure 6. Rule viewer GUI showing membership 
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Figure 7. Surface viewer for designed ANFIS 


5. RESULTS AND ANALYSIS 

The model which is analyzed under transient conditions with PI and ANFIS is a 9-MW wind farm. 
It is simulated for a period of 20 seconds and studied using phasor analysis. The simulation is carried out 
using solver ode23tb in MATLAB R2015a environment. The wind turbines have a nominal wind speed 
rating of 10m/s. The simulation is performed with applied speed of 8m/s on all the three turbines initially. 
Then, a change in wind is applied to turbine 1 at t=2sec and the speed changes from 8m/sec to 11m/sec. 
Similarly, wind speed is rammed to 11m/s for turbine 2 and turbine 3 at 4 sec and 6 sec respectively. 
The Figures (8-13) give the comparative analysis of active power, reactive power, bus voltage, rotor speed 
and electromagnetic torque with proposed and PI controller. 


N 





TURBINE3 
TURBINE 2 
à. 





¥(14,1.47) 





Active Power (MVV) 
ho 
Reactive Power (Mvar) 


rx... 





p 


TURBINET 








“0 2 4 6 8 1 12 14 «416 488 2% 0 2 4 6 8 0 12 4 16 18 2% 


Time (seconds) 


Figure 8. Comparison of output active power of 3 
wind turbines with PI and ANFIS based controller 
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Figure 9. Comparison of output reactive power of 3 
wind turbines with PI and ANFIS based controller 
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Figure 10. Comparison of bus voltage with PI and 
ANFIS based controller 
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Figure 12. Comparsion of rotor speed of 3 wind 
turbines with PI and ANFIS based controller 


ISSN: 2089-4856 o 


P (0.233.657) 


N 


2 R (5.213,1 864) 


Pitch Angle (degrees) 





L 1 1 1 L 
8 10 12 14 16 18 20 
Time (seconds) 





=) 
N 
a 
> 


Figure 11. Comparison of pitch angle with PI and 
ANFIS based controller 
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Figure 13. Comparison of electromagnetic torque of 
3 turbines with PI and ANFIS based controller 


It can be seen from Figure 8 and 9 that the peak overshoot and settling time of active and reactive 
power are reduced with ANFIS based controller in contrast to PI controller. Alternatively, the comparative 
analysis (quantitative in nature) of improved performance with the proposed controller is also given in 
Table 4. In other words, the electrical performance of the system is enhanced. Figure 10 shows the voltage 
at bus with PI and ANFIS based controller. It can be seen that voltage with the proposed controller shows 


minor improvement. 


It can be seen that for maximum wind change on turbines, the maximum pitch angle change is 3.657 
degrees and has been shown as in Figure 11. Using the developed controller, the pitch angle change is 
calculated to be 1.380 degrees. Therefore, for the same wind change, the mechanical movement of turbine 
blades is less with the proposed controller in comparison to PI controller. This improves the mechanical 


performance of the overall system. 


Table 4. Comparison of active power of wind farm 


ACTIVE POWER 
TURBINE 1: 
Peak overshoot (MW) Settling time (seconds) 
PI 3.237 12.7 
ANFIS 3.016 11.92 
TURBINE 2: 
Peak overshoot (MW) Settling time (seconds) 
PI 3.239 12.47 
ANFIS 3.011 12.18 
TURBINE 3: 
Peak overshoot (MW) Settling time (seconds) 
PI 3.227 13.72 
ANFIS 3.004 12.54 
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Figures 12 and 13 show the dynamic behaviour of rotor speed and electromagnetic torque with 
respect to time. It can be again seen that the rotor speed and electromagnetic torque are settling in minimum 
time with the use of developed controller in contrast to PI controller. The rotor speed and electromagnetic 
torque is termed as mechanical parameters of the system. 

Figure 14 shows the transient performance of rotor speed and electromagnetic torque. 
The investigation of these characteristics are necessary for analysing the performance of any 
electromechanical based system [20]. It is to be noted that the number of oscillations are reduced with 
the developed controller in contrast to PI controller. It is also to be noted that this Figure gives 
the performance of turbine 1 only. The similar behaviour is expected for turbines 2 and 3. 

At steady state, it can be seen from Figure 8 and 9 that the turbines output active power and reactive 
power settle at an operating point X and Y given as (14,3) and (14,1.47) respectively. Therefore, when 
the turbines attain steady state at fully loaded condition, they generate the maximum rated active power of 
3-MW. Further, it can be observed from Figure 12 and 13 that when the system attains the steady state 
operating point A, the ordinates are (14, 1.0047) and (14,-0.91) respectively. And these ordinates are 
replicated in Figure 13 by emitting the time axis. Therefore, in Figure 13, now the operating point A is 
having the ordinates (1.0047,-0.91). As this characteristics gives the clear understanding of oscillations. 

The detailed non-linear model of single input and single output wind power is converted into linear 
model using linearization tool of MATLAB. Therefore, the pole-zero plot of this linearized model is shown 
in Figure 15. It can be seen that all the poles lie on the left hand side of the imaginary axis and therefore 
the system is stable proving the validity of the work done. 
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Figure 14. Torque speed characteristics for turbine 1 Figure 15. Pole-zero plot for PI and 
with PI and ANFIS based controller ANHIS controller 


6. CONCLUSION 

In this paper, the transient performance (electrical and mechanical) of a 9-MW IG wind farm has 
been improved by implementing ANFIS based controller in contrast to conventional PI controller. An ANFIS 
based controller has been designed and implemented to improve system performance in contrast to 
a conventional PI controller. The study has been carried out in MATLAB/SIMULINK R2015a environment 
with solver ode23tb through phasor analysis. This improvement has been shown with the help of various 
waveforms like active power, reactive power, voltage, torque and speed. It is concluded that the proposed 
controller reduces peak overshoot as well as settling time. The stability shown by pole-zero plot, obtained 
by linearizing the detailed non-linear model proves the validity of the work done. 
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